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Disentangling the circumnuclear environs of Centaurus A: 
Gaseous Spiral Arms in a Giant Elliptical Galaxy 
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ABSTRACT 

We report the existence of spiral arms in the recently formed gaseous and 
dusty disk of the closest giant elliptical, NGC 5128 (Centaurus A), using high 
resolution 12 CO(2-l) observations of the central 3' (3 kpc) obtained with the 
Submillimeter Array (SMA). This provides evidence that spiral-like features can 
develop within ellipticals if enough cold gas exists. We elucidate the distribution 
and kinematics of the molecular gas in this region with a resolution of 4'.'4 x If 9 
(80 pc x 40 pc). The spiral arms extend from the circumnuclear gas at a radius of 
200 pc to at least 1 kiloparsec. The general properties of the arms are similar to 
those in spiral galaxies: they are trailing, their width is ~ 500 ± 200 pc, and the 
pitch angle is ~20°. From independent estimates of the time when the Hi-rich 
galaxy merger occurred, we infer that the formation of spiral arms occured on a 
time scale of less than ~10 8 yr. The formation of spiral arms increases the gas 
density and thus the star formation efficiency in the early stages of the formation 
of a disk. 

Subject headings: galaxies: elliptical and lenticulars, cD — galaxies: individual 
(NGC 5128) — galaxies: structure — galaxies: ISM 
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Introduction 



Some ellipticals show rotating disk-like features of gas and dust, expected to form as 
a result of gas accretion from the intergalactic medium, cannibalization of other galaxies, 
or re-accretion of material originally ex pelled during the major merger event that resulted 



2011; 


Young et al. 


2011 


, e.g.)- 



in the formation of the elliptic al itself (IHenkel fc Wiklindl 119971 ; lYoungl |2002| ; iDavis et al. 



Sage fc Gallettal (119931 ) proposed an evolutionary sequence 



of disk formation in dust lane ellipticals via accretion, from irregularly distributed gas to 
more evolved and settled disks. However, our understanding of the properties and evolution 
of such relatively settled disks is poorly understood yet, partly because most observations 
of the gas in elliptical galaxies are hampered by poor angular resolution and/or sensitivity. 
In particular, it is not known if these newly created disks in elliptical galaxies are uniform 
once settled, or whether deviations from axisymmetry such as prominent spiral structures 
can develop, affecting the subsequent star formation patterns, and the overall evolution of 
the disk itself. 



At a distance of D ~ 3.8 Mpc ([Harris et al.l 120101 ) (where 1" ro ughly corresponds to 



18 pc), NGC 5128 (hereafter Cen A) is the closest giant elliptical (Ivan den Bergh 



Israel 11998 



Harris et al. 



1990 



2009|). 



iarris et al .1120121 ) with a prominent dust lane along its minor axis ( llsrael 



1976 



1998 



20121) fsee Figured! left). This dust lane contains 



as traced by Ha (Nicholson et al. 19921). mid-IR (IQuillen et al. 20061). H I ( van Gorkom et al. 



Struve et al. 



2010f ). and CO (buillen et al.l Il992l ; [R^dbeck et al.l fl993t fespada et al. 



arge amou nts of gas and dust, 



The H I emission still shows unsettled gas at r>6 kpc, including tail/arm like struc- 
tures, and it is esti mated that the merger event occurred within a timescale of only 0.3 Gyr 
( jStruve et al.ll2010l ). The gas in molecular phase and the dusty component inside the inner 
few kiloparsecs are well settled within a warped disk of about 7' with an inclination of i ~ 
70° ( IQuillen et al.ll2006l ). although its detailed distribution at scales of a few hundreds of 
parsec has hitherto been poorly understood. 

The most popular model to reprod uce the large scale di s k structure is that of an 



homogenous warped and thin disk (e.g. 


Nicholson et a 


. 1992 


Kainulainen et al. 


2009; 


Struve et al. 


2010; 


Quillen et al. '< 


2010). 



1992| IQuillen et al.l Il992l . 12006 



20101 ). According to these models, 



the disk crosses the line of sight at two galactocentric radii. One of the most remarkable fea- 
tures in favor of this model is that in projection it can reproduce the pa rallelogram structur e 
seen in the IR, having a major side of ~3' (3 kpc) along a PA = 120° (jQuillen et al.l 120061 ) . 



However, the inner 1 k il opars ec of Cen A (up to a few 100 pc from the nucleus) is 
very complex. lEspada et al.l ( 120091 . hereafter Paper I) report 12 CO(2-l) SMA observations 
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of its nuclear region in a single pointing, with a resolution of 6"0 x 2"4 (100 x 40 pc), 
revealing the detailed distribution and kinematics of the molecular gas in the central inner 
kpc, including a compact circumnuclear disk of molecular gas surrounding the AGN, as well 
as outer molecular gas partly associated with the parallelogram structure. While the warped 
disk models reproduce the observed features (distribution, kinematics, dust lane appearance) 
at a large scale, Paper I shows that the molecular gas distribution in the inner kpc is not 
sufficiently in agreement with these models. A weak bi-symmetric potential was proposed in 
Paper I to explain the deviations with respect to this model. 

We present new observations using the Submillimeter Array (SMA0; IHo et alJliooi l. 
that allow us to shed more light into the complex few inner kiloparsecs. The resolution of 
the mosaic we present here is a factor of 45 higher tha n existing 12 CO (2-1) maps covering 
the dust lane, with a resolution of 23" (or ~ 380 pc) fjRydbeck et al.lll993[ ). Our mosaic 
covers an area three times larger than previously published interferometric CO maps (Paper 
I) to fully cover the entire parallelogram structure. 



2. SMA 12 CO(2-l) Observations and Data Reduction 

Cen A was observed at 1.3 mm using the SMA with 7 antennas on April 22 and 26 
2008. The field of view is characterized by a Half Power Beam Width (HPBW) of the 
primary beam of an SMA antenna of 52" (0.9 kpc). We observe a 5 pointing mosaic to 
cover 156 x 52"(or 2.6 x 0.9 kpc). The digital correlator was configured with 6144 channels 
(2 GHz bandwidth), resulting in a velocity resolution of about 0.5 km s" 1 . The receivers 
were tuned to the redshifted 12 CO(2-l) (v TCSt = 230.538 GHz) emission line in the upper 
sideband (USB), using Vlsr = 550 km s" 1 . Note that velocities are expressed throughout 
this paper with respect to the LSR using the radio convention. We used R.A. = 13 h 25 m 27!6 
and Dec. = -43°01'08"8 (J 2000) as our ph ase center (AGN position R.A. = 13 h 25 m 27!615 



and Dec. = — 43°01 / 08"805, iMa et al.lll998l ). The maximum elevation of the source at the 
SMA site is ~ 27°, forcing us to observe only under very stable atmospheric conditions, with 
zenith opacities of typically T225 ~ 0.10. In order to have a beam shape as close to circular 
as possible we used a compact configuration with longer N-S baselines. Minimum and 
maximum projected baselines were 6 m and 108 m. The maximum angular scale observable 
with the shortest baselines is 25", which considerably limits the amount of missing flux to 



lr The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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emission with angular scales larger than this. 

The two tracks were reduced independently. The editing and calibration of the data 
were done with the SMA-adapted MIR software!. 3C273 and 3C279 were used for passband 
calibration. An initial gain calibration was performed using J1337-129, which is at an angular 
distance of 30° from the target. We confirm that the continuum emission toward Cen A 
was found to be unresolved (Paper I). The gain calibration for the five pointings was then 
refined using the averaged line- free channels in Cen A itself (central pointing). Callisto 
and Ganymede were used as absolute flux calibrators. Overall, we estimate the absolute 
flux uncertainties on the order of 10%, from the comparison of both datasets. Finally we 
combined both tracks. 



An image of the 12 CO(2-l) emission was produced using MIRIAD (ISault et al.lll995l ). 
A careful subtraction of the continuum was done using line-free channels with the task 
UVLIN. The data were cleaned using MOSMEM with uniform weighting. The synthesized beam 
is 4"4x l'/9 (80 x40 pc) with a major axis P.A. = 25.3°. The task MOMENT was used to calculate 
the 12 CO(2-l) integrated flux density distribution {2a clipping) and the intensity- weighted 
velocity field distribution (clipped at 3a). The rms noise level of the integrated intensity 
map is 6 Jy beam -1 km s -1 . Figure |2] shows the SMA 12 CO(2-l) integrated intensity map 
and velocity field along the inner 3' (~3 kpc). 



3. Spiral arms in a Giant Elliptical Galaxy 

From the 12 CO(2-l) emission distribution in Figure [21 we confirm that the molecular gas 
is preferentially located along two filamentary structures resembling spiral arm like features 
to the SE and NW of the circumnuclear gas (galactocentric radius < 200 pc and position 
angle P.A. = 155°, Paper I). The larger field of view of our observations demonstrates that 
these filamentary structures extend at least l'.O, and they are curved toward the NE and 
SW, respectively, as projected distance from the nucleus increases. The approaching and 
receding sides of the disk are on the E and W, respectively, and the velocity range spans 
about 600 km s _1 , from 200 to 800 km s" 1 . 



The near side of the disk at r < 1.3 kpc is in the South ( iQuillen et al.ll2010l ). From the 
morphology of the bisymmetric spiral arm like features and the velocity field in Figure HI 
we infer that they are trailing (convex side advances). The projected angular widths of the 



2 MIR is a software package to reduce SMA data based on the package originally developed by Nick 
Scoville at Caltech. See |http : //cf a- www . haxvaxd . ed.u/~cqi/mir cook . html| 



-5 - 



spiral arm like features are ~10". Adopting an inclination of 70° (jQuillen et al.ll2010l ). their 
widths are 500 ± 200 pc. 

The s hape of galactic spiral arms is usually logarithmic in nature and independent of 



scale (e.g. ISeigar fc James! Il998l ). One of the best geometric measures to represent spiral 
arms is the pitch angle, 0, defined as the angle between the line tangent to a circle and 
the line tangent to a logarithmic spiral at a specified radius. We present in Figure [3] the 
deprojected CO(2-l) distribution, as well as a best fit representation of a logarithmic spiral 
pattern with pitch angle = 20°. We used a set of models of bisymmetric spiral arms with 
an inclination of 70°. We varied the pitch angles in steps of 10° from 10 to 60°, and found 
that <£> = 20° fitt ed best. This pitch angle is within the expected range for late type spirals 
(jKennicuttlll98ll ) for a maximum rotation velocity for Cen A of V ~ 300 km s _1 . 



In Figure H] we s how a comparison b etween the CO(2-l) and mid-IR 24/xm emission 
obtained with Spitzer (jQuillen et al.ll2006l ). with resolutions of about 6", which complements 
Figure [1] with its mid-IR 8/im emission. In general the molecular gas distribution is similar 
to that of the dust emission. 

As mentioned previously, a warped and thin disk model is usually used to reproduce the 
observed properti es at large scales (>1 kp c), especially the 3' parallelogram structure seen 
in dust emission (buillen et alJbood . lioioh . However, the morphology of CO (2-1) emission 
resembles better spiral arms rather than a parallelogram structure. This rules out the single 
contribution of a homogeneous disk, otherwise we would also see filaments to the NE and 
SW of the circumnuclear gas, as the warped and thin disk model predicts that they are of 
comparable brightn ess as the observed NW and SE filaments (see for example Figure 4 in 
Quillen et al.ll2006l ). This cannot be explained by flux loss due to missing spacings. From 
the model, these filaments should have sizes less than 10", similar compactness than the 
filaments to the NE an SW, and thus they would have been detected. 

Therefore, we interpret the two CO (2-1) filamentary structures as spiral arms. Their 
main properties are similar to those found in disk galaxies: they are trailing, the linear widths 
are a few hundred parsecs, and the pitch angle is typical of spiral galaxies. This provides 
evidence that spiral arms can develop within a giant elliptical if enough cold gas exists. 
However, note that the molecular gas properties in the spiral arms might be considerably 
different to those in spiral galaxies, as it coexists within a triaxial potential and the stellar 
surface density changes more abruptly with radius in the case of a giant elliptical. 
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Discussion 



A consequence of the compression produced by these spiral arms is that it is expected to 
trigger s tar formation (SF) o n the leading edge. This is consistent with abundant SF traced 
by Paa (jMarconi et al.ll2000l ) for example, coincident at least with the southern molecular 
spiral arm (nearest part and thus likely the least obscured). The spiral features seen in the 
12 CO emission is also associated with the distribution of certain tracers obser ved in poorer 



angula r /spectral resolution maps, such as the SCUBA 450^xm emission map in lLeeuw et al. 



(120021 ). as well as the pure rotational lin e of molecular hydro gen H 2 (J=2-0) S(0) (28.22 /xm) 
emission observed with Spitzer/IRS by iQuillen et al.l ( 120081 ). The mole cular hydrogen tran - 
sition H 2 (J=2-0) S(0) indicates the presence of gas with T ~ 200 K (IQuillen et al.l 120081 ). 
which is likely tracing photodissociation regions associated with abundant star formation. 

A small perturbation could have triggered this spirality, such as a non-axisymmetric 
weak potential (Pa per I) or a minor merger after the disk was relatively well settled (e.g. 
Purcell et al.l 1201 if ). From the H I structure and kinematics we infer that the formation of 
spiral arms took place in less than 0.3 Gyr (jStruve et al.ll2010l ). which is likely the most 
accurate measure of time at our disposal since H I is one of the most sensitiv e compo- 
nents of interaction. A lthough evidence suggests that spiral arms are transient (jSellwood 



2011 



Fovle et a 



(jD'Onghia et al. 



20111) . recent simulations claim that spiral structures may be long lived 



20121 ). 3-D self-gravitating models for late type spirals suggest that the 



pitch angle we find here would correspond to a long-lasting feature , because the density 
response tends to avoid larger pitch angles (jPerez-Villegas et al.l 120121 ). Similar simulations 
in the deep potentials of giant ellipticals would be needed to investigate the response of the 
gas. 

The study of the properties of recently accreted molecular gas in the deep potential 
wells of elliptical galaxies is a powerful tool in disentangling how galaxies form and evolve. 
Recent smoothed-particl e hydrodynamic (SPH ) cosmo logical simulations in the A Cold Dark 
Matter ( CDM) scenario (iKomatsu et al.ll2011l) are able to reproduce most of the properties 
of disks (lAgertz et al.l 120 111 ; iGuedes et al.l 120 111 ; iDomenech- Moral et al.l 120121 ) at the current 
epoch, including their fine structure morphology, the observed angular momentum, the Tully- 
Fisher relation, and the SFR to gas surface density ratio (Kennicutt-Schmidt). However, 
to understand galaxy evolution it is essential to anchor these simulations with the actual 
properties of nascent disk galaxies. The properties of disks, in particular the existence of 
spiral arms, in these systems which are in the early stages following significant accretion 
can be compared with numerical simulations. Not only can we compare the observed gas 
distribution and kinematics of recently formed disks to models, but also the effect of SF and 
Active Galactic Nuclei (AGN) feedback in the disk evolution. Moreover, with the advent 
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of ALMA, it is now possible to validate theories of disk formation as they enter a deep 
gravitational potential, performing simil ar studies in objects alo ng the evolutionary sequence 
of accreted gas in dust lane ellipticals ( iSage fc Gallettal Il993l ). from systems whose gas is 
still irregularly distributed to more evolved and settled disks such as in Cen A. 



We thank the referee for useful comments that improved the focus of this Letter. This 
research has made use of NASA's Astrophysics Data System Bibliographic Services, and 
has also made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. This research was partially supported 
by a Marie Curie International Fellowship within the 6 th European Community Framework 
Programme. 

Facilities: SMA 



REFERENCES 

Agertz, O., Teyssier, R., & Moore, B. 2011, MNRAS, 410, 1391 
Davis, T. A., Alatalo, K., Sarzi, M., et al. 2011, MNRAS, 417, 882 

Domenech-Moral, M., Martinez-Serrano, F. J., Dommguez-Tenreiro, R., & Serna, A. 2012, 
MNRAS, 2503 

D'Onghia, E., Vogelsberger, M., & Hernquist, L. 2012, ArXiv e-prints 
Espada, D., et al. 2009, ApJ, 695, 116 
-. 2010, ApJ, 720, 666 

Foyle, K., Rix, H.-W., Dobbs, C. L., Leroy, A. K., & Walter, F. 2011, ApJ, 735, 101 
Guedes, J., Callegari, S., Madau, P., & Mayer, L. 2011, ApJ, 742, 76 

Harris, G. L. H., Gomez, M., Harris, W. E., Johnston, K., Kazemzadeh, F., Kerzendorf, W., 
Geisler, D., & Woodley, K. A. 2012, AJ, 143, 84 

Harris, G. L. H., Rejkuba, M., & Harris, W. E. 2010, PASA, 27, 457 

Henkel, C, & Wiklind, T. 1997, Space Sci. Rev., 81, 1 

Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 



-8 - 

Israel, F. P. 1998, A&A Rev., 8, 237 
Kainulainen, J. T., et al. 2009, A&A, 502, L5 
Kennicutt, Jr., R. C. 1981, AJ, 86, 1847 
Komatsu, E., et al. 2011, ApJS, 192, 18 

Leeuw, L. L., Hawarden, T. G., Matthews, H. E., Robson, E. I., & Eckart, A. 2002, ApJ, 
565, 131 

Ma, C, et al. 1998, AJ, 116, 516 

Marconi, A., Schreier, E. J., Koekemoer, A., Capetti, A., Axon, D., Macchetto, D., & Caon, 
N. 2000, ApJ, 528, 276 

Nicholson, R. A., Bland-Hawthorn, J., & Taylor, K. 1992, ApJ, 387, 503 

Perez- Villegas, A., Pichardo, B., Moreno, E., Peimbert, A., & Velazquez, H. M. 2012, ApJ, 
745, L14 

Purcell, C. W., Bullock, J. S., Tollerud, E. J., Rocha, M., & Chakrabarti, S. 2011, Nature, 
477, 301 

Quillen, A. C, Brookes, M. H., Keene, J., Stern, D., Lawrence, C. R., & Werner, M. W. 
2006, ApJ, 645, 1092 

Quillen, A. C, de Zeeuw, P. T., Phinney, E. S., & Phillips, T. G. 1992, ApJ, 391, 121 
Quillen, A. C, Neumayer, N., Oosterloo, T., & Espada, D. 2010, PASA, 27, 396 
Quillen, A. C, et al. 2008, MNRAS, 384, 1469 

Rydbeck, G., Wiklind, T., Cameron, M., Wild, W., Eckart, A., Genzel, R., & Rothermel, H. 
1993, A&A, 270, L13 

Sage, L. J., & Galletta, G. 1993, ApJ, 419, 544 

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in ASP Conf. Ser. 77: Astronomical 
Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, & J. J. E. 
Hayes, 433 

Seigar, M. S., & James, P. A. 1998, MNRAS, 299, 685 
Sellwood, J. A. 2011, MNRAS, 410, 1637 



-9 - 

Struve, C, Oosterloo, T. A., Morganti, R., & Saripalli, L. 2010, A&A, 515, A67 
van den Bergh, S. 1976, ApJ, 208, 673 

van Gorkom, J. H., van der Hulst, J. M., Haschick, A. D., & Tubbs, A. D. 1990, AJ, 99, 
1781 

Young, L. M. 2002, AJ, 124, 788 

Young, L. M., Bureau, M., Davis, T. A., et al. 2011, MNRAS, 414, 940 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 10 - 




< 



ug 

1 





Fig. 1.— ie/tj Optical image of Centaurus A, showing its prominent dust lane 
(ESO/IDA/Danish 1.5 m/R. Gendler, J.-E. Ovaldsen & S. Guisard, ESO). Right) The molec- 
ular gas as traced by our SM A CO (2-1) observat ion (green), PAH and dust emission at 8 
jum observed by Spitzer (red, iQuillen et al.l l2006h . and the Chandra X-ray observations of 
the jet (blue, NASA/CXC/M. Karovska et al.). Note that the spiral and the 8/xm emission 
in this panel are within the optical dust lane visible in the left panel. The most prominent 
features of the 12 CO(2-l) emission are consistent with those of the 8 /j,m emission, but the 
CO shows the spiral arms. There is a CO absorpti on feature toward t he unresolved compact 
continuum component located toward the AGN ( jEspada et al.ll2010l ). Note that the jet is 
nearly perpendicular to the circumnuclear molecular gas at the base of the spiral arms, but 
not to the outer component. 
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Fig. 2. — Top) 12 CO(2-l) integrated intensity map. Color scale ranges from -20 to 
103 Jy beam -1 km s _1 . The rms of the CO image is 6 Jy beam -1 km s" 1 . Bottom) 
12 CO(2-l) (intensity weighted) velocity field map. Contours are placed every 50 km s _1 , 
from 200 — 800 km s _1 . The color scale ranges from 400 km s _1 up to 700 km s _1 . The size 
of the synthesized beam is shown in the lower left corner (dark blue filled ellipse). 
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Fig. 3. — The deprojected CO(2-l) distribution and a logarithmic spiral pattern with pitch 
angle <fi = 20° (black dashed line). 
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Fig. 4. — Comparison of Spitzer 24/xm emission (contours) and CO(2-l) emission (grey 
scale) . 



